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Biomolecules are often very long with a definite chirality. DNA, xanthan and poly-7- 
benzyl-glutamate (PBLG) can all form columnar crystalline phases. The chirality, however, 
competes with the tendency for crystalline order. For chiral polymers, there are two sorts 
of chirality: the first describes the usual cholesteric-like twist of the local director around 
a pitch axis, while the second favors the rotation of the local bond-orientational order and 
leads to a braiding of the polymers along an average direction. In the former case chirality 
can be manifested in a tilt grain boundary phase (TGB) analogous to the Renn-Lubensky 
phase of smectic-A liquid crystals. In the latter case we are led to a new "moire" state with 
twisted bond order. In the moire state polymers are simultaneously entangled, crystalline, 
and aligned, on average, in a common direction. In this case the polymer trajectories in 
the plane perpendicular to their average direction are described by iterated moire maps of 
remarkable complexity, reminiscent of dynamical systems. 
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1. Introduction and Summary 

It is well known that large molecules play a central role in the structure and function 
of the cell 0]. In particular, DNA, large polypeptides such as poly-7-benzyl-glutamate, 
and polysacharrides such as xanthan are long polymers with a definite and consistent 
chirality. DNA, with a chain length on the order of centimeters, must be packed into 
regions with length scales on the order of microns, a scale much smaller than the average 
polymer end to end distance in dilute solution. The packing of these molecules both in 
vivo and in vitro is of great interest. In the absence of specialized cellular structures, it is 
plausible that liquid crystalline phases could facilitate in this packing . It is known that 
bacterial plasmids can form both nematic and cholesteric liquid crystalline mesophases || . 
In addition, dinoflagellate chromosomes J3J and bacterial nucleoids |5[] exhibit cholesteric 
phases, while sperm heads and bacteriophages [[7| exhibit hexagonal columnar phases. 
Outside the cell, many mesophases arise for long and short chiral biomolecules: columnar 
phases P,pyT0[| , cholesteric phases P, |IT| , |T2|1 , nematic phases |T3| , p^| , p!5[] and recently even 
blue phases have been observed. The possibility of new phases arising in long chiral 
molecules is intriguing. 

In the hexagonal columnar phase of chiral liquid crystals, the crystalline close-packing 
of the molecules competes with the tendency for the molecules to twist macroscopically 
around each other. As in the twist grain boundary phase of smectic-A liquid crystals 
||16|1 , chirality can enter the crystal through the proliferation of screw dislocations. In close 
analogy with the physics of type-II superconductors, the screw dislocations enter when their 
energy per unit length is smaller than the energy gain from introducing chirality. When 
the chiral couplings are small, screw dislocations are excluded and a perfect equilibrium 
crystalline phase persists, in the same way that the Meissner phase expels an external 
magnetic field below the lower critical field H c \. 

In this paper we elaborate and extend the results presented in |T^,^8|. As before we 
neglect structure along the polymer backbones and consider the allowed chiral couplings in 
the hexagonal columnar phase. We find, in addition to the usual cholesteric-like term which 
favors rotation of the local polymer direction, a new coupling which favors the rotation 
of the crystalline bond order along the polymer axes. This term leads to a novel phase 
in which the polymers are braided and in which their trajectories can be described by 
iterated moire maps. Modulated chiral phases have been discussed in the context of two- 
dimensional films |T9| , pO| , |21[] of chiral tilted molecules. Our new phase is a three-dimensional 
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chiral modulated structure. Although the moire textures we find have some similarity with 
the blue phases of cholesteric liquid crystals, proposed earlier for chiral polymers [|12[ , they 
differ in that they also incorporate close-packed crystalline order. 

We propose two new liquid crystalline phases. The first is the direct analogue of the 
Renn-Lubensky twist grain boundary phase |16|] in which a sequence of polymer crystal 
regions are separated from each other by tilt grain boundaries which effect a finite rotation 
of the average polymer direction. The director (i.e. the polymer tangent) undergoes a 
cholesteric-like rotation around a pitch axis which is perpendicular to the average polymer 
direction as well as to the tilt grain boundaries. We depart here from the notation of Renn 
and Lubensky in calling these walls "tilt boundaries" rather than "twist boundaries" as 
in []16[]. We use the nomenclature "helical grain boundary" for the honeycomb networks 
of screw dislocations in the moire state discussed now. The second set of structures are 
phases characterized by a sequence of equally spaced parallel helical grain boundaries 
which cause rotations of the local bond-order. These twist grain boundaries rotate the 
bond order about an axis parallel to the average polymer direction and perpendicular to 
the dislocation walls. This new phase leads to complex polymer trajectories which are 
highly entangled and reminiscent of chaotic dynamical systems. One or more of these 
phases should appear in sufficiently chiral biopolymers when concentrated in a isotropic 
solvent. Additionally, we expect that discotic liquid crystals, which also form columnar 
phases, can exhibit these phases as well if the molecules are chiral. It may also be possible 
to induce the phases we discuss by imposing appropriately twisted boundary conditions 
on weakly chiral or achiral samples. 

In section 2 we formulate the theory of chiral polymer crystals, introducing the suc- 
cessive degrees of order which separate an isotropic polymer melt from a directed polymer 
crystal. We show how the sequence of phases progresses through spontaneous symmetry 
breaking and the non-zero expectations of order parameters. Additionally, we emphasize 
the analogy of crystals to superconductors |22[] and the concomitant analogy of rotational 
invariance with gauge-invariance. 

In section 3 we calculate the line energies of various types of dislocations using con- 
tinuum elastic theory. We argue that DNA as well as discotic liquid crystals should be in 
the type-II regime, thus allowing mixed phases with a proliferation of defects. 

In section 4 we analyze the tilt grain boundary (TGB) phase of the chiral polymers. 
Using the continuum elastic theory discussed in section 3, we analyze the effect of a wall 
of screw dislocations on the tilt and twist fields (nematic director and bond-order). We 
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estimate the lower critical chiral coupling 7 C for the TGB phase in terms of the Landau 
parameters. In addition we show in a continuum dislocation density approach that the 
TGB phase arises naturally as a low energy configuration of screw dislocations. Finally 
we discuss the x-ray structure function expected for a macroscopic sample. 

In section 5 we analyze the new moire phase. We calculate the lower critical field 7^ 
for this phase to exist as well as the effect on the tilt and twist fields of a single twist 
wall composed of a honeycomb network of screw dislocations. This phase is also shown, 
in a continuum "Debye-Hiickel" approach, to arise as a low energy configuration of screw 
dislocations. The braided structure of the moire phase contains parallel regions of double- 



twist [^3 1 as found in the blue phases of chiral liquid crystals. Unlike the blue phases, 
however, the order in most cross sections perpendicular to the average polymer direction 
resembles a perfect triangular lattice. 

Finally, in section 6 we propose a microscopic structure for the moire state. We model 
the polymer trajectories as lines woven by a sequence of moire maps. We discuss the 
preferred rotation angles for which moire patterns appear as well as the scaling statistics 
of the polymer trajectories. Upon projecting the polymers onto a plane perpendicular to 
their average directions, we study the scaling of the distance between neighboring paths, 
leading to Lyapunov exponents. In the Appendix we prove that the special rotation angles 
that produce moire maps are irrational fractions of 2tt and thus the moire state is an 
incommensurate structure: the projected polymer paths never repeat. 



2. Rotational Invariance and Free Energies 

We first derive the elastic theory of the polymer crystal, starting from phases of higher 



symmetry. Columnar crystals are solid-like in two directions and liquid-like in one [24]; we 
assume that the individual monomers in neighboring columns are out of registry, unlike 
conventional solids which are crystalline in three directions. At high temperatures or in 
the limit of extreme dilution the polymers will be isotropically entangled. We proceed by 
breaking successive rotational and translational symmetries until a hexagonal columnar 
phase is reached. All these phases need not exist for a given material; one or more of the 
intermediate phases could be skipped via a direct first order transition. 



2.1. The Polymer Nematic Phase 



When the temperature T is decreased or the density increased, an isotropic polymer 
melt can transform into a nematic phase |2f| . Orientational fluctuations are then described 



by the usual Frank free energy density: 



T^ = \ (^i(V-n) 2 + K 2 [n- (Vxn)] 2 + K 3 [nx(Vxn)] 2 } 

(2.1) 



ifi(V ± -5n) 2 + K 2 {V ± x8n) 2 + K z {d z 8n)' 



where the {K^} are the splay, twist and bend elastic constants. In the nematic phase 
we can take n w z + Sn where 8n is a vector in the plane perpendicular to the average 
director. Here, and in the following, vectors which lie in the plane perpendicular to the 
average direction of the director will be denoted as A, while full three dimensional vectors 
will be represented by bold face A. We will henceforth take the average nematic direction 
to be z. Because the polymers are long, K\ is much larger than it would be in a comparable 
short-chain liquid crystal [p6,27|, diverging when the polymers become infinitely long. In 



addition, non-linear effects are known to modify the behavior of K 2 and K% and the 
bulk compression modulus, leading to their anomalous dependence on the polymer lengths 
2^,^9[ . If the system were chiral, terms which violate parity transformation would be 



allowed. In a nematic phase the director has a discrete n — > — n symmetry which must be 
preserved in the free energy. Since under parity V — > — V, the lowest order chiral term 
is n-(Vxn). Note that since every term must have even powers of n we can chose to 
take n — > — n (vector) orn-m (pseudovector) under parity. For definiteness we take n 
to be a vector, though this has no effect on the analysis or results. We write the chiral 
contribution to the nematic free energy as 

jr* = _ 7n .(Vxn) « - 7 V ± x5n, (2.2) 

where the two dimensional pseudoscalar cross product is axb = e^aibj. One ground state 
configuration of T-a + T*^ is a cholesteric state with pitch qo = 7/K2, although more exotic 
blue phases are also possible |25 |. Nonlinearities associated with rotational invariance are 



known to modify the relation between 7 and go IH3 



2.2. Nematic-hexatic (N+6) Phase with Chirality 

At still lower temperatures or higher concentrations a new phase may occur which 



includes hexatic order in the plane perpendicular to the nematic director n |3(J . Although 
the present experimental evidence for the hexatic order in polymer nematics is sketchy, 
such phases seem highly likely in columnar systems, in analogy with the hexatic order 
possible for vortex lines in high temperature superconductors ||31|| . There is a complex 
order parameter, i/jq which characterizes the local bond order. Its behavior is determined 
by a Landau theory 

T % = i/n|n-V^ 6 | 2 + ^/i 2 |V^ 6 | 2 + ^e| 2 + ^el 4 (2.3) 

with r = a(T — Tq). In the ordered phase, for T < Tq, (i/jq) ^ and we may write 
^6 = IV'ele 6 * 6 ' 6 . Sufficiently below the transition we may neglect the fluctuations in the 
magnitude of the order parameter and employ a Landau theory for 9q. Upon suppressing 
an additive constant, Tq becomes 

Tq = l -K\(n.Ve,f + l -Ki [(W 6 ) 2 - (n-W 6 ) 2 ] « \K\{d z 9,f + \k^9qY (2.4) 

where k\ = 36 1 ipe\ 2 (h\ + hz) and = 36 1 -i/>6 1 2 ^2 are hexatic stiffnesses parallel and 
perpendicular to the nematic director and \if)Q \ 2 = a(Tg —T)/4u. The symmetry n — > — n 
discussed above for pure nematics must also hold in this phase. However, since 9q is 
measured around n, upon taking n — > — n we must also have 9^ — ► —9q. If, for instance, 
we had originally measured a change in 9q with respect to the right-hand rule using n, we 
would now, see the change go in the opposite direction, with respect to — n. Thus we may 
have any term quadratic in powers of 9q or n. Under parity, 9q will change sign if we take 
n to transform as a vector, for the same reason it changes sign under nematic inversion. 
Note that 9$ — ► — 9q is equivalent to i/jq <-> ipg. We are now able to write a new chiral term 
for the "N+6" phase. Namely 



T\ = -i^n-tyfiWe - ^V^ 6 *) = - 7 'n-W 6 « -j'd z 9 6 (2.5) 

where 7' = IT" | V^s I • I n a nematic ground state with n = z everywhere (i.e. 7 ~ 0), this 
term will favor a ground state with 9q twisting with a pitch qo = 7' / 'k\ along the ^-axis. 



A similar coupling has been considered before in hexatic smectic-S phases P3] . If both 7 



and 7' are nonzero and qo and go are incommensurate, one might expect that a chiral N+6 
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polymer melt would resemble an incommensurate smectic [33, 32] . These considerations 
also apply to short chain nematogens with the same symmetries: if a chiral nematic (or 
discotic) had the analogue of an N+6 phase two distinct chiral couplings would again be 
allowed. Table 1 gives a list of important couplings with their symmetry properties. 

There are, in addition to the chiral coupling between n and 9q, non-chiral couplings 



between n and 9q p0| , p4| . Based on rotational, nematic and parity invariance, the lowest 
order terms are 

J^ne = C(n-V06)[n-(Vxn)] + C"V0 6 -(Vxn) « Cd z 9 6 V±x8n + C'V ± 9 e xd z 8n (2.6) 

where C = C + C' and C = C' . While the final two terms differ only by a total derivative, 
they are, in principle, different in the presence of topological defects in n and 9q . 

2.3. Hexagonal Columnar Phases with Chirality 

At still lower temperatures T < T xta i (or higher concentrations) the polymers can 
crystallize in the xy-plane. To insure rotational invariance, and, in analogy with treatments 
of crystals made up of point particles |35[ , we build up the crystal out of a superposition of 



plane density waves with wavevectors G a . In the case of a triangular lattice we could take, 
for example, G Q = |G| [cos(7ro;/3), sin(7ra/3), 0] where a = 1,...,6, |G| = 4-7r/(\/3ao) and 
ao is the lattice constant of the crystal. The {G a } are the six smallest reciprocal lattice 
vectors of a triangular lattice and G a ■ z = 0. The areal polymer density in a plane 
perpendicular to the average direction may be expanded as 

6 

p « po + p a (r) exp{-zG a -r} (2.7) 

a=l 

where higher order reciprocal lattice vectors could also be included. 

Each plane wave is modulated by a spatially varying magnitude and phase p a (r) = 
|p a (r)|exp{zG Q! -'u(r)}, where u is a two-dimensional displacement field. Under a rotation 
by 9, the plane waves in ( |2.7|) change. In particular r — > r + 9 x r, and so G ■ r — ► 
G r + r-(Gx0). Thus rotations lead to a position dependent change in p a (r), namely 
Pa —> Pa exp{zG- (0xr)}. Likewise, under a global rotation about the a;-axis or y-axis by 
an angle 9 X or 9 y respectively, n^n + 0xn^z + 9 y ~x — 9 x y, i.e. Sn — > Sn + 9 y ~x — 9 x y. 
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Similarly, under a global rotation about the z-axis by 9 Z , 9§ — > 9q + 9 z . To insure rotational 
p2| , |35[| , derivatives of p Q (r) must be accompanied by the fields 6>6 and 5n . 



mvanance 



J~ xtal — 











G a • 


V x p Q 



V x p Q - z/9 6 (G Q xz) p c 



2 B 
+ 2 



V x p Q - i9 6 (G a xz) p c 



f2J 



q 6 1 G a +G / 3+G 7 =0 

+ ^" \d z p a - i (G a -Sn) p a \ 2 + -j\Po\ 2 > + c ^2 PotP(iP 1 + 0(p A a ) 

Crystalline order arises in this Landau expansion for sufficiently small b ~ (T — T xta i) so 
that ( p a ) 7^ . Due to the third order term, this transition will, in general, be first order, 
and T xta i represents a limit of met ast ability. The free energy is the sum of all the terms 
discussed earlier, namely 



F = Jd 3 X {F n + F 6 + ^ n6 + -^xtal} 



(2.9) 



The free energy ( |2.9| ) is analogous to the Ginzburg-Landau theory of a superconduc- 
tor. Here, rotational invariance dictates the couplings between Goldstone modes, Sn and 
6*6, and derivatives of p a in the same way that gauge invariance dictates the coupling be- 
tween derivatives of the BCS order parameter if) and the electromagnetic field A. The 
contribution (jF n + JF 6 + JF n6 ) mimics the field energy (VxA) 2 . The two chiral couplings 
are analogous to the coupling between the magnetic field B = VxA and the external field 
H ||36|| . In our case there are two distinct "magnetic fields" H represented by the chiral 
couplings 7 and 7'. 

Upon setting p a = \p\ exp{zG Q -u] and integrating out Sn and 6q we find that these 
fields become locked to various derivatives of the displacement fields: 



Srii « d z Ui 
9q ~ ~^ e ijdi u j 

The resulting free energy to lowest order in u and its derivatives is now 



(2.10) 



F = Jd 3 x \pu% + ^ul + K^dlm) 2 - 7 V x x5n - y'd z 9 6 } (2.11) 

where Uij = ^(diUj + djUi). Here and throughout Roman indices indicate directions only 
in the xy-plane while Greek indices indicate all three coordinate directions. In terms of 
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the parameters in (pUD, p = ^\G\ 4 \p\ 2 (A + B) and A = ||G| 4 |p| 2 (A-S). In contrast with 
polymer nematics, in the crystal the elastic constants do not diverge when to the nonlin- 
earities associated with rotational invariance |32] are taken into account. The simplified 
free energy ( |2.11| ) is similar to the Ginzburg-Landau theory in the London limit. 

The response of superconductors to an external magnetic field is determined by k, the 
ratio of the penetration length to the coherence length. As discussed in the next section, 
the parameter analogous to k in hexagonal columnar crystals is typically much greater than 
1, showing that these materials behave as Type II rather than Type I superconductors in 
their response to chirality. In other words the equilibrium ground states may contain 
defects. 



3. Isolated Dislocations and Their Energies 

When dislocations are introduced into the free energy ( |2.11| ) u(r) is no longer single 
valued. To account for this, we introduce a new variable w-yi which is equal to d^Ui away 



from the defects [37j. The free energy becomes 



F = j d 3 xp (^ WlJ ^ Wjl ^j + ^(wa) 2 + ^-(d z w zi ) 2 - 'ye ij d i w zj - 7' 'd z (^€ijWij), (3.1) 

where 6q = ^eijWij and 5rii = w z i. 

Dislocations are restricted so that the Burger's vector, b must lie in the xy plane 
and t-(nx6) = 0, where t is the unit tangent point along the dislocation line. The 



latter constraint eliminates dislocations which add a row of polymer ends |pl| , which we 
neglect. We introduce the density tensor a 7 i(r) = f dtdbt 7 bip(t,b,r), where p(t, b, r) 
is the volume density of dislocations at the point r with Burger's vector b pointing in 
the t direction. Since the dislocations do not end, V • t = 0, and d 1 a 1 i = 0. The 
constraint = t-(nx6) = n-(frxt) becomes, assuming n||z, b x t y = b y t x or in the case 
of many defects a xy = a yx . In terms of w^i this constraint reads diW z i = d z Wjj, i.e. 
V±-5n = —d z Wjj = —d z (5p/p) (dp is the fluctuation in the areal density of polymers p in a 
constant-z cross section) which creates the sort of long-range interactions that are central 
to theories of directed polymers |]38| , p9| , |27[| . 



Following reference [57], we can relate w 1 i to the density of dislocations inside a small 
area V as follows: 



ds 
qy as 



dui 
ds 

dui_ 

\JjtXj ry 



^ J t n -dS 6^ 2 (t n xr) 



(3.2) 



dS ^ e^jduWyi = - I dS^a 



T Jr 
where dT is the boundary of V and dS is a directed surface element. It follows that 



(3.3) 



In Fourier space (|3.3| ) can be solved for w^i and we have 



q) = _. 6 7M ,g M ^(q) + . g ^. (q)) 
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(3.4) 



where it is convenient to decompose the arbitrary function ^(^l) into longitudinal and 
transverse parts, 



tpi = i 



(3.5) 



so that Vi_ • ip(r) = —u and Vj_ x ^(r) = —n. The field ip describes the equilibrium dis- 
placement of the polymers in the presence of the dislocation density ct 7 i(r). The extremal 
equations which result from minimizing ( |2.11|) with respect to variations in -0 imply that 



[(2v + \)ql + K 3 q 4 z ]q ± -<iP 



2fx + X- K 3 ql 



-eimqzqiqjOtmj 



(3.6) 



and 

2 4n - 7 A* ~ ^3<Zk /i + A 

lnq ± + K 3 q z \q_ L XTp = — q z q s qitsjeimOt m j H -^—q s qiqpe s je jpu a vi . (3.7) 

The three terms on the rights hand sides are sources representing the three different types of 
dislocation. An edge dislocation lying in the xy-plane provides a source for (|3~^). The first 
source contribution in ( |3.7| ) is the only non-vanishing term for a straight screw dislocation, 
while the second is non-vanishing for edge dislocations which lie parallel to the z-axis. 
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Following [20], we can find the effect of dislocations on the the nematic twist Vi_>5n and 
bond twist d z Q%. Let u p = e pa jW a j. Then, ( [3.4| ) leads to V-u> = — Tr[a]. The definitions 



9« 



2 e ij W ij 



and 5rii 



w, 



also imply V-u; 



d y Sn x 



d x 8n v + 2d z 6 G . Thus 



2d z 6 6 - V±x5n = -Tr[a}. 



(3.J 



Since the diagonal components of represent the density of screw dislocations we have 
related the quantities we are interested in to the macroscopic defect density. Upon referring 
to ( P-lOj) , we see that ( |3.8|) quantifies the non-commutivity of derivatives of u. By choosing 
different dislocation complexions we can make one or both of d z 9§ and V±xSn non-zero 
and thus take advantage of the chiral couplings in (|2 . 1 1 ) . 

We can now solve for the strain field around a screw dislocation. A screw dislocation 
in a polymer crystal is shown in fig. [j], where we have chosen boundary conditions such that 
5n — > at z = ±oo. For a screw dislocation along the S-y X fiix8(y)S(z)b or, in 

Fourier space, a 7 i = 27vb5 lx 5i x 5(q x ). Immediately from ( [3^4]) we see that w xx = w 
Since q± ■ ip = we also have w 



■I'lj 



0. 



yy 



w 



~y 



w 



0. It follows that 
2m5{q x )bq 



yx 



and 



w 7 



2m6(q x )bq y 



1 - 



1 + 



m \ + K 3 qj 

ql((i-K 3 q 
fiq\ + K 3 q 



4 

z 



(3.9) 



(3.10) 



From these results we can readily determine the bond twist and nematic twist fields 
generated by a single screw dislocation with the boundary conditions as above. The bond 
twist is 

dz e 6 - 1 *~ - K * h 



2 9 Z Wy X 



4 



JilyV+Izz) 



(27T) 



2 H n 2 



16v^F 



A~ 1/2 | 



-5/2 



where 



The nematic twist is 



e 4X \y\ 



X= 4 ^ 



2 9y + -^3? 



2A 



(3.11) 




(3.12) 



9yW zx 



V±x8n = d x w zy 

= bS(y)S(z) + 2d z e 6 



lib 

T 



dq y dq z q y 
(2tt)2 f^ + i^ 4 ' 



(3.13) 
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Note that ( |3.11| ) and (|3.13| ) are related as in (|3.8| ). 

Boundary conditions play a crucial role in determining the integrated value of the 
twists. The key issue is whether the boundaries at large \z\ go to infinity before or after 
the boundaries at large \y\. It is enlightening to consider, for instance, a vortex disclination 
in a two-dimensional XY model. While a symmetric configuration with the spins pointing 
out radially from the core may have the lowest elastic energy, boundary conditions can 
change this configuration. In an L x L region of the xy-plane we could impose boundary 
conditions such that the spins are normal to the surfaces at x = ±L. All the winding of 
the order parameter at the boundaries will occur along the lines at y = ±L. In the case of 
a dislocation in the columnar crystal, either 9q or 5n "wind" , depending on the boundary 
conditions. The boundary condition 5n — > as z — > ±00 implies that there is no net twist 
of the nematic director. In the above example and in fig. |T] we have enforced this constraint 
on surfaces at z = ±L for all values of y. This is equivalent to taking the large y cutoff to 
infinity first, i.e. q y — > 0. Hence in this case, the integral over a constant x plane of d z 6§ is 



dydz d z 9 6 (r) 



lim lim d z 6e(q L ) 

q z ^0 q y ^0 



b 

2' 



(3.14) 



while the similar integral of V±xSn is 



dydz Vj_X(5>n(r) = hni lim Vi_x5n(q) = 0. 

q z ^Qq y —>0 



(3.15) 



The boundary conditions and the order of limits will be of importance when adding up 
the effects of a collection of defects, as we shall see in the following sections. 

Upon inserting the expressions for uiji into Q3.1| ), we find the screw dislocation energy 
per unit length 



1 f d 3 q 
Lj (2tt) 3 



A* 1 



— \w 



yx I 



+ — \q z w z 



2(2tt) 5 



dq y dq z 



1 / V25 

H 1= arctan - 

y/2 \l-5 2 



> 



(3.16) 



where 



2 \ 1/4 



(3.17) 
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and £j_ and £ z are the short-distance y and z cutoffs, respectively. In the two extreme 
limits, the elastic free energy per unit length becomes 



fs 



oo 

(3.18) 





so the deciding factor is the dimensionless parameter S. In analogy with superconductors, 
8 2 plays the role of k, the ratio of the penetration depth to the coherence length. Indeed, 
suppose for simplicity that the two cutoffs £j_ and £ z are comparable ~ £ z = £) so that 
5 2 ~ y^K s //j/^. From inspection of ( |2.11j ), we see that A = ^K s /fi is the length scale over 



which bend deformations heal, so y K% / fx plays the role of the London penetration depth 
in a superconducting analogy [23]. Thus S 2 is the ratio of the healing length of director 
fluctuations over the healing length of density fluctuations (£ ~ ao, the lattice constant) 
which is precisely the form of k, in superconductors. When 6 <C 1 the polymer crystal is 
type-I and chirality will be excluded until the crystal breaks down completely forming a 
chiral liquid. However, when 5 ^> 1 type-II behavior occurs and chirality can creep into 
the crystal through the proliferation of defects. We must determine the magnitude of 5 to 
decide which sort of behavior we expect. 

We first consider chiral polymers, such as DNA. We take £j_ to be the average spacing 
ao between the polymer strands and £ z to be the on the order of the polymer diameter, or 
alternatively the base pair spacing in DNA. Then K3 = k B TL p /£j_ and fx ~ eo/£j_> where 
L p = Kb/k B T is the persistence length (Kb is the polymer bending stiffness) and eo is the 
polymer interaction energy per unit length. With the experimental parameters in [ JI] , we 
have £ 2 w 5A, £_l = 30A, L p = 550A and e « 1.5/c b TA _1 leading to 5 w 5. Thus we 
expect type-II behavior. 

For dense chiral discotics, we take £j_ and £ z to be the respective sizes of the disk- 
shaped molecules. In this case we expect that K 3 ~ Uq/^ z and fx ~ ^0/(^^2) where 
Uq is a characteristic interaction energy. If the interactions were entirely steric we would 
expect that Uq ~ k^T. Using (|3.17| ) we then find that 5 = £j_/£ 2 - In discotics, as the 
name implies, ^± ^> ^ z and so we expect that they will also be type-II columnar crystals. 
Note that in the limit of rod like-molecules (as in nematics) £j_ <C £ z and we expect type-I 
behavior in the hypothetical situation of nematic molecules forming a hexagonal columnar 
phase. 
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Straightforward estimates, along the lines taken in [52] , show that the screw disloca- 
tion core energy per unit length, when the magnitude of the Burgers vector b = ao and 
5 > 1, is 

E c = cK 3 ^ 4 a 3 /2 (3.19) 

where c is a numerical constant of order unity. Thus E c is the same order of magnitude 
as the 5 — > oo limits of the elastic energy displayed in ( |3.18|) . 

Considering now an edge dislocation with a Burgers vector of magnitude b and parallel 
to the z axis a 7 j = 5 lz 5i x b5(x)5(y), and the free energy per unit length is 

/z _ edge = .J^L / dqxdqy 4 + Vi + + = 3£ ln (R) (3 . 20) 



2(2tt) 2 J ™ ^ q 8 16tt U 

which, unlike the screw dislocation energy, will diverge as the logarithm of the system size. 
The core energy is E c = c' lib 2 in this case. The logarithmic divergence arises because a 
vertical edge defect acts on each xy-layer as a normal edge dislocation in a two-dimensional 
crystal. 

Unlike infinitely long polymers, discotics can also have edge dislocations lying in the 
xy-plane. For an edge dislocation running along the y axis, we have = b5 iy 5i x 5(x)5(z) 
and find an elastic free energy per unit length 

_ (2/* + X)K 3 b 2 f J g 2 z 
/xy " edge " 2(2tt)2 J dqxClqz (2fi + \)ql + K 3 qi { J 

which is finite, and can be reduced to the complicated expression ( |3.16[ ) by replacing \x 
with (2fjL + A) P[. 



4. The Tilt Grain Boundary Phase 

4-1- Strains, Energies and Displacements 

The twist grain boundary state was first proposed as the analogue of the Abrikosov 
flux lattice in chiral smectic-A liquid crystals []TB| and later discovered to exist in nature 
P5[j . In analogy with this state, and, as suggested by Kleman [pH| , the tilt grain boundary 
phase of polymers is a crystalline version of a polymer cholesteric, where the cholesteric 
twist arises from a sequence of low-angle tilt grain boundaries. 

We first consider a single grain boundary. The wall is made up of screw dislocations, 
pointing along the x direction, stacked up in the xz-plane. An example with three screw 
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dislocations is shown in fig. |^. To determine the twists and strain, we superpose the strains 
from each individual screw dislocation. With screw dislocations at z = 0, ±d, ±2gL . ., we 
have 

oo 

[V_Lx5n] TGB = ^ [V±x5n(x, y, ^ - nd)]i (4.1) 



n= — oo 



where the notation [. . .]i refers to the result ( ^.13| ) for a single screw dislocation. With the 
help of the Poisson summation formula ( |4.1|) can be rewritten as mixed Fourier transform 



m =-oo ' ? (42) 

6„ <>*r» m ,iA (2ir) 2 bX 



imz/d m 2 e -\{2irym'\y\/d' 



d K »> ^ 2d 3 

m m 

where ( |3.13| ) has been used in the second line. Since d is small compared to y and z in 
the far field limit, we approximate the second sum by the first term, the remaining terms 
being exponentially suppressed: 



[V±x5n} TGB « b6(y) 2^ 6 \ z ~ nd ) ^ — cos ( —J- J e 



^ — cos 



-(27r) 2 A|y|/d 2 

(4.3) 

^ p -(2^) 2 A|y|/d 2 

Similar manipulations lead to 

(2tt) 2 6A _ ( 2tvz\ _ (27r) 2 A|y|/rf2 



[d z e 6 ] TG B « -^r 1 - cos ( ) e ~ l J M/ > ( 4 - 4 ) 

and again ( |3.8|) is satisfied, since now the density of dislocations is precisely a xx = 
bS(y) J2n^( z ~ n d) ~ (b/d)5(y). In this superposition of many dislocation lines, stacked 
along the z axis, we implicitly took the limit in which the boundary at large \z\ goes to 
infinity first: It is the q z = term in the Poisson summation formula which puts the Dirac 
5- function in the expression for V±x5n and not into d z 6§. 

Upon inserting these results into ( |2.11|) , using our earlier results for screw dislocations, 
and neglecting the interactions between the lines, we find that the threshold chiral coupling 
above which screw dislocations will penetrate is 7 C = /screw /b where /screw includes both 
elastic and core energies. For 7 > 7 C it becomes energetically favorable for screw dislo- 
cations to flood into the crystal and form grain boundaries until their density is limited 
by the repulsive interactions between screw dislocations, fig. |3| shows a schematic phase 
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diagram for the polymer crystal. The transition to a tilt grain boundary phase occurs 
along the 7 axis when 7' is small. In the following section we will calculate 7^, the critical 
chiral coupling for the bond-order chirality. 

The tilt grain boundary texture is also suggested by a continuum elastic "Debye- 
Hiickel" approach similar to that used in |35| and | 31j| . This continuum theory requires 



length scales not only large compared to the lattice constants, but also large compared to 
the dislocation spacing. We require that the dislocation density a^j lead to displacements 
with finite elastic energies. Suppose there is a configuration which only depends, say, 
on the y-coordinate. In this case it is straightforward to show, using (|3.4j ) - (|3.7|) that 
the only configurations with non-divergent elastic energies are those for which a xx (q) — > 
const, as q — > and all other components of a^j vanish in the same limit. We consider 
only configurations with no edge dislocations parallel to the polymers, as they cost a 
logarithmically divergent energy per unit length. Upon taking a v j oc (27r) 2 5(q x )5(q z )5 ux 5j x 
we see that the only no n- vanishing component of w^j is w zx = ia xx /q y which, upon 
substitution into ( p.l|) gives T = iK 2 qoq y w 2X | q= o = — 1®- XX . Note that the quadratic 
contribution to the elastic energy in (|3.1| ) vanishes in this continuum limit. We now add 
the core energies of the dislocations, leading to the free energy T bh in the continuum 
Debye-Hiickel approximation: 

■^dh = E ijk iaij a k i - 7«xx (4.5) 
where Eijki is a positive definite matrix of line energies of dislocations lying in the xy-plane 



and ctij is the spatially averaged dislocation density tensor. As discussed in |31| the core 
energies of screw dislocations have the form 

2 / 2 

EijkiOij'Oki = E screw (aii) + E screw (oij) (4.6) 

where E screw and E' screw are core energies for screw dislocations. In the case of a single 
dislocation E' screw = 0. The minimum of (|4.5|) occurs when = 7/ (2i? screw + 2£" screw ), 
and all other components of ctij vanish. This is precisely the continuum version of the 
microscopic tilt grain boundary lattice. In the same approximation, using the relation 
between w zx and a xx , we find 5n x = w zx = - [7/ (2i? screw + 2£' / screw )]6(y), where 6(y) is 
the Heaviside step function, so that across the grain boundary 5n x jumps in proportion to 
the chiral coupling 7. Since Burgers vectors are quantized in units of the lattice vectors, 
our proposed periodic array of screw dislocations is the closest allowed approximation to 
a uniform density. 
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4-2. Structure Function 

A perfect hexagonal close packed polymer crystal aligned along the z-axis will have 
six primary 5-function Bragg peaks in the q x -q y plane at a radius 4% / (\^3ao) , where ao is 
the lattice constant of the polymer crystal. 

A tilt grain boundary, being made of parallel screw dislocations, must also be parallel 
to an allowed Burgers vector. As a result, if two crystalline regions are joined together by 
a tilt grain boundary, they will be rotated with respect to each other around a common 
reciprocal lattice vector, which we take to be G = (47r/\/3ao) [0,1,0]. Thus, in Fourier 
space, the six spots will rotate around an axis that passes through one pair of diametric 
points on the original hexagon. If the illumination volume in a diffraction experiment 
contains a number of grain boundaries, this rotation will continue, laying out the Bragg 
peaks along two circles lying in the q x -q z plane of radius 2n/ao at q y = ±27r/(v / 3ao)- If 
the angle of rotation is rational there will be discrete Bragg spots around the circle, while 
if the rotation angle is an irrational fraction of 2tt the spots will form a continuous circle, 
fig. |] shows a hypothetical structure function for a TGB state with pitch axis parallel to 
y with a rotation angle of 2tt/7. A perfectly periodic array of twist grain boundaries with 
spacing d! will produce a finely spaced set of additional peaks along the q y axis centered on 
47r/\/3ao at intervals of 2it/d' (not shown in fig. If the tilt grain boundaries are more 
irregularly spaced they will nevertheless limit the range of translational correlations along 
y and all Bragg peaks will be broadened in the q y direction with a width on the order of 
27r/d', where d' is a translational correlation length. 

We expect that due to long range interactions between the grain boundaries the angles 
of rotation will lock in at rational fractions of 2ir, in analogy with the Renn-Lubensky twist 



grain boundary phase |j45[ . While a detailed calculation of d and d', in terms of the Landau 
parameters is difficult, we can estimate their sizes by assuming that repulsive interactions 
between screw dislocations lead to d ~ d' . If the Burgers vector has length b = ao, the 
angle of rotation across each grain boundary is = tan _1 (ao/d), and thus the pitch is 
P = 2nd' /(f). Taking d ~ d! ^> ao we have P ~ 2ndd' /a®. Upon taking a typical cholesteric 
pitch to be on the order of P = 5000A and the interpolymer spacing to be ao = 50A, we 
find that d « d' « 200A and w 14°. 
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5. The Moire Phase: Continuum Elastic Theory 

5.1. Strains, Energies and Displacements 

We have seen that a TGB induces a finite jump in the director and only ripples in 
the bond-orientational order parameter which integrate to zero. We would like to find a 
dislocation structure which induces a finite jump in 9q, and thus exploits the chiral coupling 

i- 

An attractive possibility is a honeycomb lattice in a constant z-plane composed of 
screw dislocations on its links (see fig. ®. Such a structure generalizes for three-dimensional 



systems the grain boundaries discussed in [20]. To calculate the strain fields we must 



superpose the effect of a collection of finite-length screw dislocations with Burgers vector 
b. Consider a collinear row of segments, pointing along x, each of length d, separated 
by 2d. We can build up a hexagonal lattice by translating, rotating and superposing the 
strains or twists from this distribution of defects. 

Consider the dislocation density generated by the row discussed above, namely a 7 i = 
bSj X 8i x 8(y)8(z) J2 n l®( x ~ 3n<i + d/2) — <d(x — 3nd — d/2)]. The only nonzero component 
is, in Fourier space 

= 26sin(M/2) y e , iqxmd = b S in( qx d/2) y + 
Qx Qx 

where we have used the Poisson summation formula. When we Fourier transform back to 
real space, this becomes 



) y ™(W3) e<aTOa>/ M (52) 



a xx (x,y,z) = b5(y)5(zj , 

n 

We neglect the oscillatory contributions, and approximate the sum by the n = term, 

a xx w h}8{y)8{z). (5.3) 

Thus the broken line of screw dislocation segments is equivalent at long wavelengths to 
a solid line with 1/3 the Burgers vector if we are at distances larger than d. In this 
approximation, we can replace the honeycomb lattice with a triangular lattice made by 
extending all the edges to meet at the center of each hexagon, provided we divide the 
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superposition by 3 to restore the correct dislocation density. Proceeding as in the TGB 
case, we have 



[V±x5n] 



moire 

b I TT 
6AV AA 



me 



-y/nm/\A\z\ 



sm 



Tim 7T 

AA |2|+ 4 



imrjjfj_ I A 



(5.4) 



3A V AA 



7T . 7T 

AA N + 4 



cos [27T77j ■ rj_ / A] 



with ffj = z x Sj , where e*j are the unit lattice vectors of the triangular lattice of polymers 
and A = dV3/2 is the spacing between parallel Bragg planes in that lattice. 
Now, 8 z 9q has the 5-function: 



b 3 1 
[0*0 6 ] moir6 « -g^ $^(z) + -[V ± xSn 



moire 



(5.5) 



a.yy = b5(z)/(d\/?>), so Tr[at] 



To confirm ( |3.8j ), note that for a honeycomb lattice, ct X; 
2bS(z)/dVS = bS(z)/A. 

The honeycomb lattice allows us to exploit the r y' term in ( |2.11| ). The integral of 
8 z 9q over space from a stack of honeycomb lattices separated by a distance d! along the 
z-axis is —b/ (2 Ad') per unit volume. The total length per unit volume of honeycomb 
is 2d/d 2 \^3d' . Thus, assuming d,d' ^> ao so that interactions among dislocations are 
negligible, Y c = 2f screw /b. The region in which we expect the moire state is shown in fig. [| 
In fig. [B| we show a set of polymer trajectories and the dislocations leading to them. The 
moire state in fig. [5] contains an inner region of polymers which consist of a single polymer 
at the center of a bundle of six polymers twisting around it. This texture has double twist 
in the nematic field as found in the low-chirality limit of blue phases of cholesteric liquid 
crystals [p3| , |4^| . The moire state takes advantage of both double twist energies and the 
new chiral coupling 7'. The bond order field 9q and the nematic director n are linked to 
each other geometrically - since 9q is measured around n, when the director is not uniform 
we must carefully define 6q. This leads to the same sort of considerations which occur in 
He 3 — A: the wavefunction is defined in the plane perpendicular to the nuclear spin axis. 
The appropriate covariant derivative, which take the nematic "curvature" into account, is 
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(5.6) 



where f2 is the connection and its curl is unambiguously given by the Mermin-Ho relation: 

1 



(5.7) 



The geometrical connection between the moire state, the Mermin-Ho relation and the 



nearest- neighbor packing of polymers is discussed in |32[ . 

The honeycomb network of screw dislocations described above can also be understood 
in terms of a continuum approach where one looks for dislocation densities which lead to 
strains which only depend on the z-coordinate. As in the TGB phase we neglect edge 
dislocations parallel to the polymers and set a z j = 0. Upon assuming akj proportional to 
(2ir) 2 5(q x )5(q y ), we find that w z j = and 



w k i = — 



a 



yx 



yy 



-a 



xy 



k 
3 



x,y 
x,y 



(5.8) 



Upon substituting this result into ( |3.1|) , we find the Fourier transformed free energy 



2q 



2 i a y x 
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xy\ 
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2 , 2 
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xy 1 ^yx 



+ EijkiacijOtki 



(5.9) 



\piyy -\~ Oi X x\ |q=0 



In order to eliminate the terms that diverge as q z — > 0, we must take a 
edge dislocations) and a xy = a yx = a xx — a y% 



a xy [i.e. no 



0. Thus a xx = a yy are the only non-zero 



components of a^. The minimum of ( |5.9| ) occurs at a xx = a 



where we have used Q4.6|). Moreover, 9q = \ {w 



xy 



- W 



yxi 



yy 

h'/m 



7 7(8^screw + 4^ crew ) 



v + 4 J B / screw)]e(^) 

and so the bond order jumps discontinuously as we go across the helical grain boundary 
by an amount proportional to the imposed chirality. As in the TGB phase the honeycomb 
array of screw dislocations is the closest approximation to a uniform sheet distribution of 
screw dislocations in terms of quantized defects with discrete Burgers vectors. 



5.2. Structure Function 

The helical grain boundary will rotate the six primary Bragg spots of the structure 
function of a perfect crystal in Fourier space. In this case all six spots will be swept 
around a ring in the q x -q y plane of radius \q±\ = 47r/(\/3ao)- If the rotation angle is 
a rational fraction of 2n there will be a discrete set of spots, while if the angle is an 
irrational fraction of 2-k the spots will form a Bragg ring. While we might expect that 
interactions between grain boundaries would favor rational lock-in angles, we shall see 
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that the geometry of an isolated helical grain boundary favors an irrational angle. In fig. |6| 
we show a schematic structure function for the irrational angle </>3 = 2 tan _1 (A/3/21) = 9.4° 
obtained by superimposing the structure functions of 10 crystalline regions (i.e. 9 grain 
boundaries). An infinite periodic array (spacing d!) of twist grain boundaries would lead to 
a finely spaced set of Bragg peaks along the q z axis, at positions q z = 2nn/d'. If the twist 
boundaries are spaced more randomly, the translational correlations along z will be still 
limited to a range of order d' . Reflecting this fact, the Bragg spots will be broadened along 
the q z axis, with a width « 2n/d', as indicated in fig. |6|. When many incommensurate twist 
boundaries are included in an illumination volume, the diffraction rods in fig. ^ will merge 
into a continuous ring, broad along q z , but with a very narrow radial width. We have 
verified this numerically by calculating the scattering from a perfect moire state composed 
of 24 grain boundaries. This Fourier transform is axially symmetric about q z and in fig. |7| 
we show the structure in the q y -q z plane. 

From equation ( |5.5| ) we see that across a helical grain boundary the angles changes 
by (p ~ ao/2A (with b = ao the minimum Burgers vector). As in the TGB case, an exact 
calculation of d and d' would require a detailed accounting of elastic interactions between 
screw dislocations. We can again estimate their size by taking d ~ d' ^> ao- In this case 
the pitch is P = 2iv\/3dd' /ao- If a typical pitch is P = 5000A and ao = 50A then we find 
that d ^ d' = 150A and « 11°. 

6. The Moire Phase: Microscopic Model 

6.1. Iterated Moire Maps 

When one rotates one lattice with respect to another around a common point, it is 
known that at certain angles there will be a non-zero density of points coincident to both 
lattices By choosing the rotation angle across the honeycomb dislocation network to 
produce a high density of coincidence lattice sites j|6| we produce especially low strain 
energies across the boundary. The superposition of triangular polymer lattices below and 
above the boundary forms a moire pattern, fig. || shows a sequence of four especially 
simple angles for a triangular lattice. Polymers in the lower half-space (circles) must be 
connected to the closest available polymer in the upper half-space (crosses) to minimize 
bending energy. The map has a discrete translational symmetry, in the sense that any 
coincidence site could be a center of rotation. Note that an identical moire pattern could 
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have been obtained by rotating about a point of 2— or 3— fold symmetry of the coincidence 
lattice. Especially simple moire maps arise for the rotation angles 



2 tan 



v/3 



3(2n + l) 



n 



1,2,.. .. fig. |8] shows these maps for n = 1, . . . , 4, with angles <pi ~ 21. 8 C 



(6.1) 



13.2 C 



(ps ~ 9.4° and ^ 4 » 7.3 C 



It is shown in Appendix A that all such angles are irrational fractions of 2% [IS 



so 



that the structure never repeats upon iteration. Around each coincidence point there are 
n concentric rings of helical polymers. The lattice of coincidence points is also a triangular 
lattice, but with a spacing a n = ao v/l + 3(2ra+ l) 2 /2, where ao is the original lattice 
constant. The geometrical origin of such energetically preferred lock-in angles has no 
analogue in chiral smectics. The exact choice of lock-in angles and spacing between moire 
planes must be settled by detailed energetic calculations. For a square lattice an especially 
simple sequence of lock-in angles is given by 



/^square 



tan 



2n + 1 
2n(n+ 1) 



(6.2) 



leading to coincidence lattice spacings a^ quare = ao^n 2 + (n + l) 2 [18 



7 = 7 



A schematic of the helical grain boundary energy E(<p) as a function of angle for 
is shown in fig. ||. The deep cusps correspond to the special lock-in angles 



discussed above. See [[|7[ f° r a discussion of similar phenomena in conventional crystals. 
When 7' 7^ the energy has the form 



E(</>) = E{<j>) - c 7 



(6.3) 



where c is a constant. As 7 increases, the preferred minimum eventually jumps from the 
origin to one of the preferred lock-in angles. In the continuum approach of Section 4, 
we estimated a moire twist of </> ~ 11°. The n = 3 moire state nicely approximates this 
continuum result with ^3 m 9.4°. In conventional crystals, cusps in the energy landscape 
near a lock-in angle (j) n will have a contribution from the density of extra dislocations, 
proportional to — n as well as the logarithmic interaction between the dislocations, 
leading to the form Ae ex \<j> — (j) n \ In \<j> — <j> n \. In the columnar crystal, due to the softer 
interaction between dislocations, the leading term in the energy will behave as Ae oc \<j>— <j> n \ 
which is still sharp enough to produce the lock-in angles. 
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It is interesting to consider the effect of the moire mapping on nearest neighbor bonds. 
Consider a polymer crystal at z = — oo. In the xy-plane the polymers sit at the sites of a 
triangular lattice. We may draw lines connecting nearest neighbors as shown in fig. [10|(a), 
with polymers at the vertices. After allowing various arrangements of screw dislocations to 
pierce the crystal, we can look again at the polymers in the xy-plane at z = +00. Again, 
the polymers will sit at the sites of a triangular lattice, typically rotated with respect to 
the original one. The effect of intervening defects shows up if we undo the rotation but 
retain the original nearest neighbor connectivity of the polymers at z = —00. fig. pT0|(b) 
shows the effect of a single screw dislocation, which has "sheared" the polymers on either 
side with respect to the other. Note that each vertex still has the coordination number 



six it had at z = —00. As shown in fig. 11, when moving a bond, in any particular 



plaquette, we may make one of two types of moves denoted by o\ and 02 • A row of o\ 
or o"2 moves in a triangular lattice represents the effect of a left or right handed screw 
dislocation, respectively. Note that both moves are area preserving, in that the area of 
the parallelogram remains unchanged. We may represent the effect of a plane of defects 
by making a sequence of moves. The restriction of keeping a coordination number of six 
at each vertex amounts to the restriction that dislocation lines may never end and that 
they must meet in triples at 120°. fig. 0(c) shows the effect of a section of a helical grain 
boundary, near the intersection of three defect lines. 

Upon two iterations of the moire map separating three regions of polymer crystal, the 
first coincidence lattice is rotated with respect to the second coincidence lattice by precisely 
the angle of rotation <p n . Thus the composite coincidence lattice is the "coincidence lattice 
of coincidence lattices", with lattice constant a\/a$. Moire maps iterated p times lead to 
triangular composite coincidence lattices with spacing a n (a n /ao) p ~ 1 , i.e. to ever sparser 



lattices of fixed points with intricate fractal structure in between them. fig. [12] shows the 
projected polymer paths for a lock-in angle of a square lattice (with (pi = tan _1 (3/4) »s 
36.9°) iterated p = 1, . . . , 4 times. Note the intricate fractal structure which appears after 
several iterations of the map. 

6.2. Polymer Trajectories and Lyapunov Exponents 

In contrast to the TGB state polymers braided by moire maps with p>l are highly 
entangled and wander far from straight line trajectories, fig. [5| shows the polymer trajecto- 
ries for the n = 1 moire map, of a triangular lattice iterated nine times. In fig. |13] we show 
forty random polymer trajectories after 99 iterations. Although these polymers are clearly 
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influenced by their proximity to a special center of rotation, exceedingly complex trajec- 
tories are superimposed on their slow drift around this center. The center becomes less 
and less noticeable for distant polymers. Using a perfect moire state with 3721 polymers, 
we have calculated the average Fourier transform of the monomer density p(q±,q z ) of a 
single polymer, averaged over all members of the array. In fig. [14] we show contour plots 
of the intensity ( \p(q±, q z )\ 2 ) in both the q x -q y and q y -q z plane. Again we see that struc- 
ture is approximately axially symmetric about q z . If a dilute concentration of deuterated 
polymers were introduced into the mix, neutron scattering should produce these incoherent 
averages. In addition, we have studied numerically the scaling behavior of the polymer size 
as measured by the components of moment of inertia tensor in the radial and azimuthal 
directions relative to the center as a function of N, the number of iterations. We expect 
that the radius of gyration of the projected polymer, R should be a function of N and 
i?o, the radius of the starting point of the polymer trajectory. We postulate for the radial 
radius of gyration that 



Ro \ Rq r 

For large Rq or small N, we find that the radial R r and azimuthal R^ radii of gyration do 
not depend on Rq and both scale as N^, i.e. as a simple random walk in the xy-plane. 
Thus as z — > 0, A r (z) ~ z and x r + y r = |. Similarly, the azimuthal radius of gyration 
scales as 



Rq \ Rq J 

and again we find that x$ + y<f> = \ and A^[z) ~ z as z — > 0. These observations suggest 
the one parameter scaling form 

g?Wg°) = T (JL\ ( 6 . 6) 



where r\i = 1/yi, i = r, <fi, and the scaling functions Ti{z) ~ z as z — > 0. fig. |15| suggests 
that, superimposed on the small "random walk" behavior is a slow drift around the 
rotation center in the azimuthal direction. This drift represents the effect of the distant 
exceptional coincidence site. As the polymers spread around the circle at radius Rq they 
appear to stay confined to an annular region of width ~ Rq. For simplicity we take 
r\ = r\ r = rj^. We choose n so that the points of crossover to ballistic (for R^) and constant 
(for Rr) behavior collapse. We find 77 = 1 gives the best collapse. In fig. [15| we show the 
radius of gyration data for 1364 polymer trajectories around a supercoincidence site of 99 
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n = 1 moire maps. We have binned the data into 30 equal width annuli to and assigned 
an average -Ro to each bin. The figure clearly shows that ccj + yi = |. After the crossover 
points Rr(N, Rq) — > constant and R<f>(N < Rq) ~ iV, as expected. In the radial case, there 
is a second crossover point in R r (N, Ro), presumably where the polymers have made a 
complete circle around the supercoincidence site. This second crossover is evident in the 
plot only for small values of Rq. 

Not only do the polymers spread from their starting point, but they also diverge from 



their neighbors, fig. 16 shows seven polymers which start as nearest neighbors. After 99 
iterations they separate from each other, as shown by their trajectories. Consider a pseudo- 
dynamics under which the z-axis becomes time and we consider the trajectories of the 
polymers as the trajectories of two-dimensional particles. We can look at the average square 
distance between pairs of neighboring polymers and consider the scaling of the separation 
as a function of N, in the spirit of measuring a Lyapunov exponent in a real dynamical 
system. We have performed this numerical experiment with the same set of 1364 polymers 



used to study the radius of gyration. We calculate A = a/ ( \Sfi — Sfj\ 2 ) = yA^ + A|, 
where Sfi is the displacement of polymer i from its starting point, and where the average 
is over all % and j which are nearest neighbors. Again, we postulate scaling forms for the 
root mean square radial and azimuthal separations A r (N, Rq) and A ( j ) (N, R ): 



Rq \ Rq 

From our analysis we find that, again, for small N, Aj(iV, Rq) does not depend on Rq. 
Thus Bi(z) — > z as z — > and, according to our data, Xi + jji = 0.8. Thus, we will have 
the scaling form 

where = l/jji, i = r,(f). We find the best collapse of data for u r = = 0.3. The 



collapsed data is shown in fig. [T7|. The separation between neighbors, for small N is not 
quite ballistic. Presumably this is due to the fact that the polymers stay with their partners 
for some time before splitting apart as seen in the four rightmost paths in fig. |1^. 
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Appendix A. Irrationality of Lock-in Angles 

The special lock in angles that we have found for the triangular lattice are all irrational 
fractions of 2n. We will prove a more general theorem in this appendix that applies to any 
rotation which leads to a coincidence lattice with a finite lattice constant. An alternate, 
algebraic proof appears in JT3] . 

Consider two lattices with the second one rotated with respect to the first by an angle 
around a lattice point in common (so that if = the lattices coincide). At certain 
"magic angles" cp n there will form a lattice of points which are on both lattices. This is 
the coincidence lattice. Take the magic angle in question to be </>* . If the angle is not an 
angle of symmetry of the lattice (e.g., 27rn/6 for a triangular lattice) then the spacing of 
the coincidence lattice will be larger that the spacing on the original lattice by a factor 
of C- Now consider adding a third lattice, rotated again by 4> with respect to the second 
lattice (i.e. rotated by 2(f) from the first lattice and by (p from the second lattice). The first 
two lattices create a coincidence lattice and the second two lattices create a coincidence 
lattice. However, the last two lattices may be obtained by rotating the first two lattices 
by (p* . Thus the coincidence lattice of the last two will be rotated by </>* with respect 
to the coincidence lattice of the first two. Therefore the coincidence lattice of all three 
lattices will be the coincidence lattice of the two coincidence lattices and so the spacing of 
the supercoincidence lattice will be a factor of ( 2 larger than the spacing of the original 
lattice. 

This argument can now be extended to the coincidence lattice of n lattices. The 
first (n — 1) lattices form a coincidence lattice, as do the last (n — 1) lattices. These two 
coincidence lattices are rotated by (p* from each other since the last (n — 1) lattices can be 
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obtained by rotating each of the first (n — 1) lattices by 0*. Thus the coincidence lattice 
of n lattices is the coincidence lattice of two (n — 1) coincidence lattices. This means that 
every time we add a another rotated lattice, the coincidence lattice spacing grows by a 
factor of ( and so the coincidence lattice of (n + 1) lattices is a factor of ( n larger than the 
original lattice spacing. 

If the angle 4>* were a rational fraction of 2tt, 4>* = 2-7r|, then after no more than 
q iterations the lattices would start to repeat. Thus upon adding the next lattice, the 
supercoincidence lattice would remain the same, and thus the lattice spacing would not 
grow. This implies that C g_1 = 1 an d so the angle <p* must be an angle of symmetry of 
the lattice. This proves that the lock-in angles are not rational fractions of 2tv. 
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Operator 


aymmeiiy unciei 
Spatial Inversion 


oymmeny unciei 
Nematic Inversion 


v = V6> 6 


+ 


_ 


n 


— 




vn 




+ 


v(Vxn) 


+ 


+ 


n-(v-V)n 

v. (rt - X7\r\ 
V I 11 V III 

n(n-V)v 






n-(Vxn) 




+ 


(v-n) [n-(Vxn)] 


+ 


+ 



Table 1. Symmetries of a variety of operators. Those with parity (— ) are chiral. Only combi- 
nations which are invariant under nematic symmetry are allowed. 
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Figure Captions 



Fig. 1. A single screw dislocation in a polymer crystal. The dark horizontal line is the 
screw dislocation. 

Fig. 2. View of a tilt grain boundary, looking down the y-axis. For clarity we only show 
the polymer rows immediately behind and in front of the TGB. The heavy lines 
pointing in the x direction are the screw dislocations. 

Fig. 3. Phase diagram of a chiral polymer crystal. Insets are representative tilt (TGB) 
and moire grain boundaries. Shaded lines are screw dislocations. Although we 
focus here on TGB and helical moire states, more exotic screw dislocation phases, 
including those with melted dislocation arrays, could also appear. 

Fig. 4. Structure function of a tilt grain boundary phase in Fourier space. Due to the 
interruption of order by grain boundaries spaced along the y-axis (with period 
d'), the Bragg spots are broadened along the g^-axis. This schematic shows a 
TGB phase with rational rotation angle (p = 2ir/7. 

Fig. 5. The moire state. The thick tubes running in the z direction are polymers, while 
the dark lines are stacked honeycomb arrays of screw dislocations. The intersec- 
tion of these polymers with any constant z cross section away from the hexagonal 
defect arrays has the topology of a triangular lattice. 

Fig. 6. Structure function of a moire state in Fourier space. Because of the periodicity 
along the z-axis, the Bragg spots are broadened along q z . This shows a moire 
phase with irrational rotation angle </>3 ~ 9.4° with 10 crystalline regions. 

Fig. 7. Numerically computed scattering contours in the q y -q z plane for a moire state 
with 24 helical grain boundaries. We have scattered from each of 96 beads evenly 
spaced along each of the 3721 polymers. The structure function is axially sym- 
metric around the q z axis. This numerical evidence supports the schematic of 

fig-B 

Fig. 8. A single moire map for triangular lattices with n = 1,2,3,4 (cp n = 21.8°, 13.2°, 
9.4°, 7.3° respectively). The crosses can be thought of as the ends of polymers 
in a perfect crystalline region, as can the circles. They join by connecting to the 
nearest polymer in order to reduce the elastic energy. The shaded lines are the 
screw dislocations which make up a honeycomb network. 

Fig. 9. Change in energy from the ground state of an untwisted crystal as a function of 
the angle of rotation across a grain boundary, after Balluffi, Komem and Schober 



47 . We expect cusps in the energy around each lock-in angle to grow as 
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In a conventional crystal the cusps are sharper, behaving as — <j> n \ In — <j> n \. 
Nonetheless, the cusps in the energy landscape will be responsible for the lock-in 
to the special angles <f> n . 

Fig. 10. Connectivity diagram for a polymer crystal, (a) shows the unperturbed connec- 
tivity at z = — oo. (b) shows the effect of a single screw dislocation. Note that 
every vertex has coordination number six on a cross section at z = +00, preserv- 
ing the nearest neighbor connectivity of the polymers at z = —00. (c) shows a 
section of a honeycomb array of screw dislocations near an intersection of three 
dislocations. Again the coordination number of every vertex is six. 

Fig. 11. Allowed moves on a triangular lattice. Note that the moves preserve the area 
between the four vertices. These moves can be put together to model the effect 
of dislocation arrays. 

Fig. 12. The projected top view of a moire map on a square lattice with rotation angle 
tan _1 (3/4). The four boxes show the projected polymer paths after p = 1, 2, 3, 4 
iterations. 

Fig. 13. A projected top view of 40 random polymers paths resulting from the moire map 
with n = 1 iterated 99 times. There is an exceptional fixed point of all 99 maps 
at the center. 

Fig. 14. Numerically calculated average Fourier transform of the monomer density of a 
single polymer in the (a) q x -q y and (b) q y -q z planes. We expect that a dilute 
concentration of deuterated polymers would produce neutron scattering profiles 
similar to these. The plots were made from the same data set as in fig. [7[ 

Fig. 15. Radial (R r ) and azimuthal (R^) radii of gyration near a supercoincidence site. 

The polymers are binned according to Rq, the distance of their initial point 
from the center of rotation. We have rescaled the radii and N by Rq and col- 
lapsed all the data. The lower curve is log 10 [R^(N) / Rq] and the upper curve is 
log 10 [ J R2( iV )/ J R ] + 2. 

Fig. 16. Trajectories of seven nearest neighbors. The large dots represent the starting 
points of each polymer path. 

Fig. 17. Radial (A r ) and azimuthal (A^) root mean square separation of nearest neighbors 
near a supercoincidence site. The polymer groups are binned according to Rq, 
the distance of their initial point from the center of rotation. We plot the rescaled 
values of A versus log 10 [A^/i?Q- 3 ]. The lower curve is log 10 [A r (A^)/i?Q' 24 ] and the 
upper curve is log 10 [A^ ) (A^)/ J R[]- 24 ] + 2. 
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